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Abstract 

This report concerns the energy of a zero-temperature many-body system of spin \ fermions in the 
unitary limit. In a previous report (nucl-th/0705.0944) this energy was determined to be £ = 0.24 in units 
of the free gas kinetic energy, appreciably lower than most reports giving £ ~ 0.45. In our calculation the 
2-body interaction satisfied exactly the unitary limit i.e. infinite scattering length and effective range ro = 0. 
In the present report results with ro > are shown. A strong dependence on the effective range is found. 
It is for example found that an increase to ro = 1/m increases £ to ~ 0.4 close to other reports of £ in the 
unitary limit. It is concluded that because of the singular character of the unitary limit it is necessary to 
verify that the interaction actually satisfies unitarity. 

The calculations done here in a pp-ladder approximation show a resonance in the in-medium interaction 
close to (and in) the unitary limit. This was already found in the previous work. 
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1 Introduction 

P^P . The properties of a dilute fermigas with large scattering length is of considerable theoretical as well as experi- 

mental interest. 

Several numerical methods have been used to determine £ with results varying from £ ~ 0.5 to £ ~ 0.25. 
Extensive references are for example found in [TJ [2] 

Apart from the interest in the properties of a unitary gas these calculations also provide a test of many- 
body methods. The Monte Carlo calculations are at least in principle the most accurate and could provide a 
benchmark. 

It is a well-known fact that interactions with large scattering lengths are separable. This paper is a report 
on results of calculations using separable interactions determined by inverse scattering from phase-shifts of 
various large scattering lengths a s and small effective ranges including a s = oo and r = .0. The energy at 
, zero temperature is calculated from an in-medium effective interaction obtained from a pp-ladder summation 

J> ' of these separable interactions. 

Section 2 is a short summary of the method used with numerical results shown in Section 3 and some of the 
conclusions are summarised in Section 4. 
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2 Separable Interaction from inverse scattering 

The method used here to calculate a separable interatction by inverse scattering has also been used in several 
previous papers where details can be found. 3, 4j[5j[6] Only a short summary of the expressions is given here. 
The input in the calculations are phase-shifts, either experimental or otherwise defined as shown below. A rank 
one separable potential provides a sufficient and in fact precise description of the interaction in the unitary 
limit. If the phase-shift changes sign such as in the experimental 1 Sq case a rank two potential is necessary. 



(see ref [3]). In the case of a rank one attractive potential one has 

V(k,p) = -v(k)v(p) (1) 

Inverse scattering then yields (e.g. ref [31(7]) 

v 2 (k) = ^L s in8{k)\D{k 2 )\ (2) 
k 



where 



D(k 2 ) = exp 



tt J k 2 -k' 2 ^ 



(3) 



where V denotes the principal value and S(k) is the phascshift. A provides a cut-off in momentum-space. The 
effect of the cut-off will be exploited below. 
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With 5(k) = 7r/2, the unitary limit, one finds 



and vl(k) -» -oo for k -> A. (See Fig [2]). 
For A>i one finds 



(A 2 -fc 2 )^ 



(47T) 5 



A 



(4) 



(5) 



In this limit, but only in this limit , the unitary interaction will then be a (5-function in coordinate space. And 
the strength is inversely proportional to the cut-off. 

The Viow k approximation is adequate at low density for some of the cases shown below but NOT in the 
unitary limit. An effective in-medium Brueckner interaction as defined by a particle-particle ladder summation 
is used for all the results presented here. Dispersion corrections are expected to be small [6] so that the 
denominator has only kinetic energies; the effective mass m* = 1. The diagonal elements of the in- medium 
interaction is then 



G(k,P) = - 



v 2 (k) 



l + I G (k,P) 



with 



Ia(k,P) 



(2tt) 3 



(0) 



(7) 



where P is the center of mass momentum and Q the angle-averaged Pauli-operator for pp-ladders (Brueckner 
approximation) . One should note that the angle-averaging is exact here because the denominator is independent 
of P. 

It was shown in ref [4] that this can be rewritten as 



G(k,P) = - 



v 2 (k) 



I GK (k,P) 



with 



I GK (k,P) = 



(2tt) 3 



2kf 2q ^(fc^)-^ ^,, , kv 2 (k) 

V ( fc ) k 2_ k ,2 k dk 



(8) 



(9) 



tan<5(fc) 

Eqs (JS1 EI) have the advantage over eqs (JHJ [TJ) in that the integrand in eq ^ is zero for k' > 2k f because the 
factor Q(k',P) — V is then equal to zero. Consequently, there is no need to resort to a low-momentum (" Vi ow 
approximation here. Although the two sets of equations are numerically identical the latter set simplifies the 
computing greatly. This is in particular important in the unitary limit in which case the integration can be 
done analytically to give eq (fTTj) . 

With a = yj and y — the potential energy per particle PE/A is 



PE/A 



8G(a,y)y 2 dy + 



(l-a 2 )7 



4G(a,y)(l-y 2 - a 2 )ydy 



a da 



(10) 



The kinetic energy per particle, i.e. the uncorrelated fermi-gas energy is given by 

3 h 2 



E FG /A=——kj. 

10 m 1 



The total energy is expressed in these units by 



Ej A - ZE FG /A. 



In the unitary limit (S(k) = n/2) there is a simplification already used in ref [4]. The cutoff A can then be 
chosen large so that eq ([5]) is valid and with v u {k') independent of k' one can perform the ^'-integration in eq 
(J9j) analytically. After dividing by v 2 one then finds 



Iu(a,y) 



1 + y + a * log 



1 + y — a 



+ -L(l-y 2 -a 2 )log 

2y 



(i + y) 2 -' 



l-y 2 -a 2 



(11) 
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Figure 1: This figure shows the phase-shifts that were used as input for the calculations with the labels defined 
in the text below. 



and 

G(a,y) = -47r[I u (a,y)}- 1 (12) 

This then provides an analytic expression for G(a, y) and as already shown in ref [4] one finds in this case 
£ = 0.24, explicitly independent of A and of density (fermi-momentum) as is to be expected in the unitary limit. 

It is of interest that Steele [5], using the effective field theory power counting method arrives at the same 
expression for G(a,y). But in eqfTUlfor the potential energy he uses the approximation G(a,y) ~ G(0,0) and 
then arrives at £ = 4/9. 

3 Numerical Results 

The expression §§§ was used in all calculations together with the analytic expression (jlip . valid the in unitary 
limit. 

Calculations were made with phase-shifts defined as follows: 



A: The experimental 1 Sq phases. 



B: Scattering length a s 
C: Scattering length a s 
D: Scattering length a s 
E: Scattering length a s 
F: Scattering length a s 
G: Scattering length a s 



— 18.5/m and effective range ro — 2.8/m. 
— oo and effective range ro = 2.8/to. 
— oo and effective range ro = 1.0 fm. 
— oo and effective range ro = 0.5/m. 

— oo and effective range ro = O./m and the interaction given by eq. (j4|). 
• — oo and effective range ro = O./m and the interaction given by eq. (|5j). 



Fig [T] shows these phase-shifts as a function of momentum for k < A = 4/m. Notice that phases "A" and 
" B" are comparable for small momenta, that " B" and " C" overlap except for small momenta and that " F" and 
"G" are for the unitary limit, S = ^ . 

The potentials correponding to each of these sets of phase-shifts were calculated using the inverse scattering 
eqs © and ([3]) with the results shown in FigO The potential curves A — > D with effective ranges r < 1.0 
are (practically) overlapping for k < lfm~ 1 even though the scattering length varies over a large interval, 
— oo < a s < —18.5. Note however the drastic change of the potentials when r changes from 1 to 0.5 to 
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Figure 2: This figure shows the diagonal in momentum-space of the potentials corresponding to the phase-shifts 
in Fig[Q 



(curves labelled D,E and F) while a s — > — oo. Note that curve "F" being the exact solution for a separable 
potential in the unitary limit is uniquely different from all the others. 

The total energy £ in units of the fermigas-energy EpcjA is shown in Fig [3] Our curve "A" for the 1 So 
phases agrees practically exactly with the corresponding curve of ref |2J even though the latter also includes 
some ring diagrams. Curve "B" with the nl So" scattering length and effective range gives a somewhat larger 
value of £. Note that, with a s held constant at — oo, £ changes from ~ 0.4 — » 0.5 for ro = 2.8 to ~ 0.24 for 
ro = 0. The value £ = 0.24 independent of kf is indicated by the straight line. It was already obtained in ref [1] 
using eq Ijlip. Curve "F" and more so, curve "G", nearly coincide with this analytically obtained result. Curve 
"G" uses the same large A limit of v u given by eq (O as in the analytic calculation, while "F" is calculated 
with A = 6/m . The value of £ = 0.24, given by the straight line, is the unitary limit result from the analytic 
calculation. It is only in this limit that the analytic integration is possible. The numercal results represented 
by curves "F" and "G" are necessarily less accurate as discussed below. Those calculations were done only in 
order to confirm computational consistency. 

It was already pointed out in ref [4] that, in the unitary limit, Iqk = along a line in the (a, y)-plane from 
~ (0.8,0) to ~ (0.9,0.3). The associated singularity in G(a,y) implies a resonance with the in-medium phase- 
shift — > 7r/2 which appears very unusual. It may however not survive in a higher order calculation, in particular 
including spectral broadening as in Green's function calculations. The singularity in G(a, y) complicates the 
numerical integration in eq (|10[) and an interpolation method was used in ref [1] . The same method was now used 
for obtaining the "F" and "G" results for £. Instead of the analytic expression (fTTj) the numerically integrated 
eq ([9]) has now to be interpolated and this results in less accuracy resulting in the oscillations seen in Fig [3] 
especially for curve " F" . Because the analytic solution is more accurate it supersedes the " F" and " G" and an 
improvement in the numerical calculation to remedy this situation was ruled unnecessary. 



4 Conclusions 

It is , we believe, conclusively shown that a serious study of a system in the unitary limit has to be done in the 
limit itself. Any extrapolation is questionable, especially as regards the effective range. The unitary limit is (in 
a broad sense) singular and unique. This is the necessary conclusion from Fig [3] It is in fact suggested already 
by studying the phase-shifts shown in Fig [1] and even more so the potential-curves in Fig [5] There is no reason 
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Figure 3: This figure shows the energy (in units of the fermi-gas energy) for various fermi- momenta. The 
somewhat "wavy" appearance of the curves "F" and "G" (the unitary limit) is because of numerical difficulties 
described in the text. The straight line is the explicitly density-independent result from an analytic solution of 
the effective interaction G in the unitary limit. [4] 



to believe that this conclusion is unique to the methods used here. The strong dependence of £ on the effective 
range tq is an important effect and the major result of this investigation. 

Using eqs ([8]) and ([9]) the cut-off can be chosen arbitrarly large, A 3> 2k f and eq (O is then a valid 
approximation in the unitary limit. This presents a great simplification resulting in the analytic expression (|12|) 
for the effective interaction. This same result was obtained by Steele[5] using power-counting. Eqs (|6l7p are 
numerically equivalent with (|8I9[) but requires using eq ^ for the interaction resulting in a much more difficult 
integration with the upper limit equal to A where the unitary interaction Q diverges. 

As for the numerical value of £ it is not claimed to be decided by the zero-temperature ladder approximation 
used here. Dispersion (mean field) corrections wre not included as only kinetic energies were used in the 
propagators, i.e. m* = 1. This approximation led to eq ([9]) and eq which provides a great simplification 
and indeed made the calculations feasible. The dispersion correction is small for the 1 So state [6] and although 
not proven it is expected to be small also in the unitary limit. Dispersion-corrections are related not only to the 
mean field but also to the in- medium two-body corrrelations. [B] . It is in general repulsive and would therefore 
increase the value of £. It can not be ruled out that it is larger than expected. 

Pairing corrections are important. The critical temperature is relatively high (e.g. [T]) and the boson-fraction 
is high at zero temperature. 

It is generally accepted that Monte Carlo calculations is the best approach to solving the unitary problem. 
The present investigation is presented only to suggset the importance of using an excatly " unitary interaction" . 

The calculations were based on using separable interactions obtained by inverse scattering from phase-shifts. 
Earlier work has shown the near equivalence of this method with the conventional potential approach [3]. It is 
therefore expected that the conclusions obtained here are not unique to our use of these separable potentials 
but are in fact general. 

In a future publication a report with Green's function techniques allowing a finite temperature calculation 
including hh-ladders will be reported. These calculations are more extensive than the ladder summations used 
here. The effects of spectral broadening as well dispersion effects are included. 

As already pointed out above there are some similarities of these calculations with those of ref [2] although 
that work only considers r$ > 2.54/m for which there is however a semi- quant at ive agreement between the 
respective results even though a type of ring-diagrams are included in ref [5] and the interactions are also 
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different. 

After this report was first submitted to the arxiv the author became aware of a work by T. Schaefer et al [5] 
extending Steele's work. One result of their work was the same value £ = 0.24 referred to above and reported 
in ref. [4]. They also obtain a strong dependence on effective range as reported here. 

I thank Professor Schaefer for kindly bringing this publication to my attention. 
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